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The bremsstrahlung process v T — > v e e + e~y is the dominant photon-producing decay of heavy r 
neutrinos with m„ T > 2m e so that this process is instrumental for supernova 1987A constraints on 
the properties of these particles, notably on their mixing amplitude with v & . We calculate the rate 
of the bremsstrahlung process as well as the photon spectrum and the angular distribution relative 
to the spin direction of the parent neutrino. We carefully pay attention to the difference between 
Dirac and Majorana neutrinos. 

PACS numbers: 14.60.Pq, 14.60.Lm, 13.35.Hb, 97.60.Bw 
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I. INTRODUCTION 

The direct experimental limits on neutrino masses 
are so crude that astrophysical arguments remain the 
most important source of information on their possi- 
ble magnitude. The well-known cosmological limit 
m v < 40 eV improves the experimental m l , T limit of 
about 18.2 MeV || by almost six orders of magnitude, 
assuming that it is stable on cosmological time scales. 
This longevity cannot be taken for granted because the 
standard-model decay v r — > v e e + e~ is fast on cosmic 
time scales even for relatively small mixing angles be- 
tween v T and v e . Therefore, to plug this loophole on the 
cosmological m VT limit, one needs restrictive limits on the 
v e -v T mixing angle. Even then, of course, heavy r neu- 
trinos could decay fast by new interactions so that the 
cosmological mass bound is never absolute. In fact, heavy 
z/ T 's with "invisible" fast decays can fix certain problems 
of the standard cold dark matter cosmology, leading to 
the interesting class of cosmological rCDM models || . In 
any case, proving that a heavy v r requires interactions 
beyond the particle-physics standard model to escape the 
cosmological mass limit is a significant constraint on its 
possible properties. 

If v e and v T mix with each other, any v e source also 
produces a certain v r flux. Experimental limits on 
the flux of decay electrons and positrons from beam- 
dump reactor and solar neutrinos || provide 
strong limits on the mixing amplitude. The most re- 
strictive limit arises from core-collapse supernovae which 
are powerful v T sources because they produce this fla- 
vor directly rather than by an assumed mixing with v e . 
The decay positrons would linger in the galaxy for about 
10 5 years until they annihilate so that one can derive 
restrictive limits on the v T — > i> P e + e~ channel from the 



observed interstellar positron flux 0-^] . In addition, the 
absence of 7-rays associated with the neutrino burst from 
supernova (SN) 1987A provides further limits where the 
bremsstrahlung process v T — > i/ e e + e~j is the dominant 
photon-producing reaction ■ 

This literature has two important gaps. First, no de- 
tailed calculation of v T — > v e e + e~~f is available. While 
the simple estimates used in Refs. [p|-pd| are certainly ad- 
equate for a first estimate, a detailed calculation of the 
expected 7-ray spectrum from a supernova depends on 
the angular distribution and spectrum of the photons in 
a given v T — * v e e + e~^i decay. It is the purpose of our 
paper to provide this missing calculation. 

A second gap concerns the data used for the deriva- 
tion of the limits. Refs. |8]-[l0|] used 7-ray data that had 
been taken by the Gamma Ray Spectrometer on the So- 
lar Maximum Mission satellite in coincidence with the 
SN 1987A neutrino signal. Ref. jll[] used 7-ray data 
from the Pioneer Venus Orbiter spacecraft that were also 
taken in conicidence with the SN 1987A burst. How- 
ever, to constrain radiative decays from MeV-mass neu- 
trino decays it is not critical to use data even close to 
the i> e burst because the expected 7 pulse is widely dis- 
persed. The COMPTEL instrument aboard the Comp- 
ton Gamma Ray Observatory satellite looked at the SN 
1987A remnant for about 0.68 x 10 6 s in 1991 Q, much 
longer than the previous experiments. Because of the 
long COMPTEL viewing period, one could derive lim- 
its on the v — > 1/7 channel which are far superior to 
the SMM and PVO results, assuming that the neutrino 
mass exceeds about 0.1 MeV [|l2|. For smaller masses 
the 7 burst would have ended before COMPTEL looked 
at the SN so that the previous limits remain of interest. 
In the v T —> v e e + e~ scenario it is necessarily assumed 



that 



> 1 MeV so that the COMPTEL limits are 
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naturally the most restrictive. However, the authors of 
Ref. |l2j have only analysed the direct channel v T — * v x j. 
Our bremsstrahlung calculation could be used to extend 
the interpretation of the COMPTEL data to the more 
interesting v T — > f e e + e~7 case. 

Aside from the practical application in the supernova 
context, the bremsstrahlung process turns out to be an- 
other interesting example for the differences between 
massive Dirac and Majorana neutrinos. 

We proceed in Sec. II with a brief discussion of the 
bare process v T — > v e e + e~ . In Sec. Ill we calculate the 
rate, spectrum, and angular distribution of the brems- 
strahlung process v r — > is e e + e~"f. In Sec. IV we discuss 
and summarize our results. 
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FIG. 2. Phase-space factor according to Eq. (0). 



II. THE BARE PROCESS 

We begin with a discussion of the bare decay process 
i> 3 — > i>ie + e~ where we assume a neutrino mass hierarchy 
of the form 7713 > m 2 > m\. The 7773 eigenstate is taken 
to be the dominant admixture of the v r flavor eigenstate, 
and similarly m\ is the dominant mass component of the 
electron neutrino. This process is then essentially what 
we sloppily called v r — > v e e + e~ in the introduction. 

For a Dirac neutrino the relevant amplitude is the 
tree-level graph of Fig. |l|(a) which was evaluated by 
Shrock jl3| . The decay rate is found to be 
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where XJ e \ and U e % are the mixing amplitudes between v e 
and V\ or 1/3, respectively, and Gf is Fermi's constant. 
It was assumed that m\ is negligibly small. The phase- 
space factor (Fig. p|) is 
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with a = (rn e /m3) 2 . $ approaches unity for 777,3 ^ We- 
This process was also studied by Li and Wilczek 14 
with an emphasis on the difference between Dirac and 
Majorana neutrinos. Following their treatment we note 
that for Majorana neutrinos the decay rate is effectively 




(a) (b) 
FIG. 1. Decay amplitude for a heavy neutrino. 



the incoherent sum of the two Dirac graphs depicted in 
Figs. [|(a) and (b). While Li and Wilczek focussed on 
the angular distribution of the final-state charged leptons 
and did not explicitly compare the absolute rates, their 
treatment implies that the Majorana decay rate is twice 
the Dirac one, a result which we also find by an indepen- 
dent evaluation of the relevant S- matrix element. The 
phase-space factor of Eq. (||) is the same in both cases. 

If 7/3 is a Dirac neutrino the normalized positron and 
electron spectra dN/dE± = T^^ ic dT\ iaze ,/dE± arc 
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where q 2 = m\ — 2mzE±+m 2 and f(q 2 ) = (1 — m 2 e jq 2 ) 2 
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FIG. 3. Positron and electron spectrum for V3 — * u\e + e~ 
according to Eq. (^) for a Dirac neutrino with 7773 ^> m e . 
The dotted line is the average between the two spectra which 
pertains to both e + and e~ when 1/3 is a Majorana particle. 
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For a parent neutrino mass 7723 3> m e these spectra are 
shown in Fig. |[ For v 3 the role of electrons and positrons 
is interchanged. For a Majorana parent the electron and 
positron spectra are the same. Each is the average of the 
Dirac spectra of Eq. (H) — see the dotted line in Fig. ||[ 



III. BREMSSTRAHLUNG 



A. Matrix Element 



For a Dirac neutrino the amplitude for v 3 — > vie + e~"f 
is given by the Feynman graph of Fig. |l](a) with a pho- 
ton attached to either the electron or positron line. We 
calculate the squared matrix element, summed over all 
final-state spin and polarization states. However, be- 
cause relativistic Dirac neutrinos are produced primar- 
ily with negative helicities we do not average over the 
initial-state helicities — a hclicity-averaged rate would not 
pertain to a realistic source. We find 
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where C\+ = P[P' + and so forth with P{ + _ 1 = Pi, + ,_, 7 
the four momenta of v\, e + , e~ , and 7, respectively. Only 
for the parent neutrino we have P3 = P3 — 7773 S3 with 
S3 its spin four vector; we have not summed over its 
helicities. In the ^3 rest frame P3 = 7713(1,-83) with S3 
a unit vector in the spin direction of 1/3. 

If the initial state is a Dirac antineutrino we find the 
same squared matrix element with the role of e + and e~ 
interchanged. Moreover, we then have P3 = P3 + 7773 S3, 
leading to P3 = 777.3(1, +S3) in the 7/3 rest frame. A simple 
CP transformation reveals that the photon angular distri- 
bution relative to the spin-polarization vector is opposite 
compared to that of a neutrino. Put another way, neu- 
trinos and antineutrinos with the same momentum but 
opposite helicities will produce identical angular photon 
distributions relative to the momentum direction. 

For Majorana neutrinos the squared matrix element is 
the sum of Eq. (Q) with the same expression under the 
exchange of e + with e - . In this case the photon angular 
distribution in the parent rest frame is isotropic, indepen- 
dently of its helicity state. For Majorana neutrinos, then, 
one may average over the initial-state neutrino helicities 
(one may use P3 = P3 ) without loss of generality. 



B. Photon Spectrum 

In order to calculate the photon spectrum we need 
to integrate over the four-body final-state phase space, 
leaving only the dio integration undone, where cj is the 
photon energy. The shape of the photon spectrum, in- 
tegrated over all emission angles, is the same for Dirac 
neutrinos or antineutrinos as well as for Majorana neutri- 
nos, independently of their polarization state. We have 
performed the phase-space integration numerically by a 
Monte Carlo technique with the parametrization outlined 
in Appendix [a]. 

To derive a useful representation of the resulting spec- 
trum we borrow from a standard result for the classical 
spectrum of an inner bremsstrahlung process Jis) ] 
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where u is the photon energy, Tbare the rate for the bare 
process, a ~ 1/137 the fine-structure constant, and v the 
speed of the produced charged lepton. This calculation 
applies to a process like beta decay where the recoiling 
nucleus does not contribute to the bremsstrahlung rate. 
In our case we have two charged leptons so that we ex- 
pect twice this rate, for the moment ignoring interference 
effects. Of course, the 1 /w divergence represents the well- 
known infra-red bremsstrahlung behavior. 

Motivated by this classical result we represent the pho- 
ton spectrum by 



dT 



— Ag(x), 
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where g(x) is a dimensionless function which depends on 
the parent neutrino mass 7773. The dimensionless coef- 
ficient A is factored out such that 5(0) = 1, i.e. g(x) 
describes the spectral shape, not its absolute amplitude. 
The factor Tbare includes a global factor of 2 for Majo- 
rana neutrinos as discussed in section II. The dimension- 
less photon energy x = lo /w max is in the range < x < 1 , 
where 
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is the end point of the photon spectrum. 

In the upper panel of Fig. we show g{x) for several 
values of 777,3. F° r & U cases g (x) decreases monotonically 
from 1 at x = to at x = 1. This behavior suggests a 
power-law representation 



(8) 



where the power-law index p(x) is itself a slowly varying 
function of x. We show it in the lower panel of Fig. [| for 
several values of 7773. In a practical application it would 
probably be enough to use a constant, typical value for 
p. For reference we show in the upper panel as a dotted 
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FIG. 4. Upper Panel: Numerical results for the dimen- 
sionless bremsstrahlung photon spectrum g(x) as defined in 
Eq. (g) for parent neutrino masses m-j, = 1.1, 2, 10, 30, and 
100 MeV (bottom to top). The dotted line is a power law 
of the form Eq. ^ with p = 1.75. Lower Panel: Power-law 
index p{x) as defined in Eq. (^[) for the same parent masses 
(top to bottom). 

line g{x) for a fixed value p = 1.75 which would appear to 
produce a reasonably good approximation if 7713 is taken 
to be around 10 MeV. 

The dimensionless coefficient A measures the zero- 
point intercept of the spectrum. Numerical values for 
several parent masses are given in Table |. 

It is a useful check of our calculation to compare these 
numerical values with the classical result Eq. (|5|) which 
should be exact in the limit of soft photons. However, 
it does not include the interference between the radia- 
tion emitted by both final-state charged leptons. What 
we list in Table | as the classical prediction, then, is 
the incoherent sum of the expected radiation power from 
the final-state electron and positron where we have inte- 
grated over the velocity distribution of these particles as 

TABLE I. Values for the coefficient A as denned in Eq. 



1713 [MeV] A (quantum) A (classical) 

100.0 7.99 7.55 

30.0 5.58 5.16 

10.0 3.44 3.06 

2.0 0.76 0.54 

1.1 0.068 0.036 



given by the bare process in Eq. (||). This classical result 
is always smaller, indicating that the interference effect 
is constructive on average. In any case, for high parent 
masses the importance of interference disappears as one 
would have expected because then the bremsstrahlung 
photons are primarily emitted co-linear with the charged 
leptons so that the overlap of the two "bremsstrahlung 
cones" becomes small. The most important conclusion is 
that we do not seem to have lost factors of 2, ir, and so 
forth in the course of our calculation. 



C. Angular Distribution 

The photon angular distribution is isotropic for Ma- 
jorana neutrinos while it is nontrivial in the Dirac case. 
A possible angular distribution can arise only from the 
terms = (P3 — m^S^P-y in the squared matrix ele- 
ment Eq. (||). Because it is linear in C*3 7 it is clear that 
after all unobserved phase-space degrees of freedom have 
been integrated out we are left with a term which is in- 
dependent of the angle 6 between the photon momentum 
and the spin-direction, and one that is proportional to 
cos6>. (The signs are chosen such that cos = +1 cor- 
responds to the direction along the parent spin.) Thus 
for a fixed photon energy oj the most general normalized 
angular distribution is of the form 

dN 1 1 + a cos 

Jcose ~ 2 ' 

where —1 < a < +1. Again, because 9 is measured 
relative to the spin direction, we have a — > —a for an- 
tineutrinos. For Majorana neutrinos we have a = for 
all to as mentioned above. 

In Fig. I we show a as a function of the dimensionless 
photon energy x — w/o; max for several choices of parent 
neutrino mass. Contrary to what one might have ex- 
pected the angular distribution depends very sensitively 
on the photon energy. Soft photons are emitted nearly 
isotropically, with a slight bias in the spin direction. The 
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FIG. 5. Photon angular-distribution parameter a(x) as 
defined in Eq. (H) for several choices of parent neutrino mass. 
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highest-energy photons are emitted dominantly against 
the spin direction if 7713 3> m e . In the limit — ► 00 it 
is easy to show analytically that a = — 1 for x = 1. (The 
reverse for antineutrinos.) 



IV. DISCUSSION AND SUMMARY 

We have calculated the rate, photon spectrum, and 
photon angular distribution of the inner bremsstrahlung 
process v T — > v e e + e~^ which could have produced a 7 ray 
burst from SN 1987A, contrary to several observations. 
Therefore, one can derive limits on the relevant neutrino 
mixing angle even though the most relevant data, the 
1991 COMPTEL observations of SN 1987A ||, have 
not yet been used in this regard. 

After assuming a certain magnitude and spectrum of 
the primary v T and v T flux, a prediction of the brems- 
strahlung photon flux and spectrum depends on the pho- 
ton energy distribution as well as their angular distribu- 
tion relative to the neutrino momentum. 

The simplest case is that of Majorana neutrinos where 
the angular distribution in the neutrino rest frame is 
isotropic, independently of the parent polarization. The 
rough estimate of the bremsstrahlung rate that had been 
used in previous papers is dT/dio = (cx/tt) Tbare/w with 
photon energies oj up to half the parent neutrino mass. It 
was overlooked that for Majorana neutrinos T^are is twice 
that of Dirac neutrinos. We have a further factor 2 be- 
cause bremsstrahlung arises from two final-state charged 
leptons. Finally, one gains the factor A which is given in 
Table |[ For example, taking m„ T = 10 MeV one gains 
overall more than a factor of 10. However, the spec- 
trum decreases somewhat faster with uj than had been 
assumed. 

Therefore, the previous limits on the mixing angle will 
improve accordingly. They would improve much further 
by applying our results to the COMPTEL data which 
encompass a much larger viewing period than the previ- 
ously used SMM and PVO data. It would an important 
future project to analyse the COMPTEL data with re- 
gard to the v T — ► i^ e e + e~7 decay! 

Dirac neutrinos exhibit a nontrivial photon angular 
distribution if they are polarized. Because the iv T 's and 
£V's emitted from a SN have opposite polarizations, the 
anisotropy of the photon distribution does not cancel be- 
tween them. On the other hand, for a mass in the 10 MeV 
range and typical energies of around 30 MeV they are not 
particularly relativistic so that their degree of polariza- 
tion is incomplete. For smaller masses the degree of po- 
larization is stronger, but the deviation from isotropy is 
less pronounced. A detailed treatment of the Dirac case 
is probably too complicated to be worth conducting in 
view of the many uncertainties with regard to the overall 
neutrino flux and spectrum. 

We do not believe that interpreting, for example, 
the COMPTEL data requires an implementation of our 



bremsstrahlung results in every detail. (One would need 
numerical results on a finely spaced grid of parameters.) 
Rather, we think that our results should be used as a ba- 
sis for an approximate treatment where the magnitude 
of the possible errors is understood and controlled. 

Independently from the SN application, the brems- 
strahlung process features differences between Dirac and 
Majorana neutrinos which are interesting in their own 
right. 
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APPENDIX A: PHASE-SPACE INTEGRATION 

We give some details about the integration of the four- 
particle phase space in the v% — > v\e + e~~i process. Fol- 
lowing Ref. Jl(| a four-body phase-space integral can be 
transformed according to 

I — 1 \ 2—1 / 

_ _7T_ f V /A(s, s a , s b ) ^X(s a ,ml,ml) 
256 J s s a 

X ^ X ( Sb ' ^ ^ ds a ds b d COS 6 dQ a dn b , (Al) 
Sb 

with the four momenta Pi — (-E,,p$), i = 1, ... ,4, and 
the masses mf = P 2 . Further, we use 

A(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2ac - 26c. (A2) 

The four particles have been grouped into two subsystems 
with s a = (P1+P2) 2 and s b = (P3+P4) 2 while s 2 = (P r + 
P2+P3 + P4) 2 = m? ni . The angles ft Qi6 are within the CM 
frames of the two subsystems while 9 is the angle with 
which the two subsystems move in opposite directions 
relative to some chosen direction in the overall CM frame, 
i.e. in the rest frame of the decaying particle, for example 
relative to its spin. The limits of integration are 

(mi + m 2 ) 2 < s a < (y/s — ma - m^ 2 , 

(m 3 + m 4 ) 2 < s b < (Vi-yil) 2 . (A3) 

Finally we need to express the original four momenta Pi 
in terms of the new integration variables. Assuming that 
the two subsystems move in opposite directions along the 
z-axis we find in the CM frame 
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The boost-factors are 
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For the other subsystem we have 
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